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weights  can  also  be  made,  incrementing  the  weights  of  all  points  above  a  given 
threshold  in  one  of  the  coordinates.  The  weights  are  assumed  to  be  always  mono¬ 
tonlc  in  both  coordinates.  An  efficient  implementation  of  this  structure  is  pre¬ 
sented  and  two  main  applications  are  described.  The  first  is  to  the  problem  of 
optimal  assembly  of  code  for  computers  with  two  kinds  of  jump  instruction:  long 
and  short.  The  task  in  the  second  application  is  the  implementation  of  a  queuin 
discipline  based  on  the  ranks  with  respect  to  two  different  criteria.  — 


Dynamic  Monotone  Priorities  on  Planar  Sets"*^ 

(Extended  abstract) 

Michael  J.  Fischer  Michael  S.  Paterson 

Yale  University  University  of  Wam'ick 

New  Haven,  Connecticut  Coventry,  England 


Abstract 

A  monotonic  priority  set  is  a  new  data  structure 
which  supports  maximum-finding  and  deletions 
over  a  set  of  weighted  points  in  the  plane.  Global 
updates  to  the  weights  can  also  be  made,  incre¬ 
menting  the  weights  of  all  points  above  a  given 
threshold  in  one  of  the  coordinates.  The  weights 
are  assumed  to  be  always  monotonic  in  both  co¬ 
ordinates. 

An  efficient  implementation  of  this  structure 
is  presented  and  two  main  applications  are  de¬ 
scribed.  The  first  is  to  the  problem  of  optimal 
assembly  of  code  for  computers  with  two  kinds  of 
jump  instruction:  long  and  short.  The  task  in 
the  second  application  is  the  implementation  of  a 
queuing  dbcipline  based  on  the  ranks  with  respect 
to  two  different  criteria. 

1  Monotonic  Priority  Sets 

Let  P  be  a  multiset  of  points  p  =  (*^y)  in  the 
plane.  P  is  partially  ordered  by  the  product  or¬ 
dering  defined  by  (xo,yo)  <  (xi,yi)  iff  xq  <  xi 
and  yo  ^  yi-  Each  point  p  has  an  associated 
weight,  w(p).  This  real- valued  weight  function  w 
is  monotonic  in  the  sense  that  if  pi  <  P2  then 
tr(pi)  <  We  define  the  following  opera- 

^This  work  wts  supported  in  part  bj’  the  Office  of  Naval 
Research  under  Contract  Number  N00014-82*K-0164 
and  by  a  Senior  Fellowship  from  the  Science  and  En> 
gineering  Research  Council. 


tions  on  P. 

•  FIND_MAX  returns  a  point  in  P  of  maximum 
weight. 

•  DELETE(p)  removes  p  from  P. 

•  V_UPDATE(yt),  a)  adds  a  to  the  weight  of  each 
point  (i,  y)  with  y  >  yo. 

•  H_UPDATE(a:o>  a)  adds  o  to  the  weight  of  each 
point  (x,  y)  with  x  >  arp* 

We  assume  that  the  updates  preserve  the  mono- 
tonicity  of  w.  (This  is  assured  if  ever}'  o  is 
non-negative.)  A  data  structure  that  supports 
these  operations  is  called  a  monofontc  priority 
get.  In  the  next  section,  we  describe  an  im¬ 
plementation  of  a  monotonic  priority  set  on  n 
points  in  which  the  total  cost  for  m  operations 
is  0  ((n  -I-  m)  log(n  -1-  m)).  The  remaining  sections 
discuss  applications  of  this  data  structure  and  a 
complementary  -completeness  result. 

2  Efficient  Implementation  of 
Monotonic  Priority  Sets 

A  naive  implementation  of  monotonic  priority  sets 
simply  keeps  a  list  of  all  points  not  yet  deleted 
along  with  their  current  weights.  Each  operation 
is  performed  with  an  entire  scan  of  the  list. 

A  first  idea  for  gaining  efficiency  is  to  associate 
the  points  with  the  leaves  of  a  balanced  binar}’ 
tree  and  to  store  at  each  leaf  the  weight  of  the 
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corresponding  point  and  at  each  internal  node  the 
maximum  weight  from  its  subtree.  This  permits 
FIND_MAX  and  DELETE  to  be  performed  in  time 
O(logn).  Each  H_UPDATE  and  V_UPDATE  how¬ 
ever  might  affect  the  weights  of  0{n)  leaves  and 
hence  require  that  much  time  when  done  in  a 
straightforward  way. 

One  way  to  update  groups  of  leaves  at  once  b 
to  place  an  offset  value  V’j  at  internal  node  t  to  be 
regarded  as  an  increment  to  all  weights  in  subtree 
i.  Thus,  the  current  weight  at  any  node  b  its 
initial  weight  plus  the  sum  of  the  offsets  on  the 
path  to  the  root.  If  the  leaves  were  sorted  by 
their  x-coordinates,  then  an  H_UPDATE  could  be 
effected  by  altering  offsets  on  at  most  log  n  nodes. 
Similarly,  V_UPDATEs  would  be  easy  if  the  leaves 
were  sorted  by  their  y-coordinates. 

For  a  pair  of  points  (xo,yo)  and  (xi,yi)  incom¬ 
parable  in  the  product  ordering,  xo  <  xi  iff  yo  > 
yi .  Any  independent  (i.e.  pairwbe  incomparable) 
set  of  points  can  therefore  be  ordered  (from  left  to 
right)  so  that  the  x-coordinates  are  non-decreasing 
and  the  y-coordinates  are  non-increasing,  permit¬ 
ting  both  H_UPDATE  and  V_UPDATE  operations 
to  be  performed  efficiently  as  described  above. 

We  choose  for  our  independent  set  I  a  max¬ 
imal  independent  set  of  (undeleted)  points  that 
are  maximal  in  the  product  ordering.  Since  P  b  a 
multiset,  I  may  not  be  unique,  so  we  specify  fur¬ 
ther  that  among  the  possibly  several  copies  of  the 
point  (x,  y),  we  choose  one  of  largest  weight  for 
inclusion  in  /.  By  the  monotonicity  condition,  I 
contains  some  point  with  maximum  weight  among 
all  of  the  undeleted  points  in  P.  Thus,  FIND_MAX 
and  the  two  updates  can  all  be  performed  on  thb 
structure  in  time  O(logn). 

However  following  a  DELETE  of  some  element 
from  the  set  I,  new  elements  may  need  to  be 
brought  into  1  to  restore  the  maximality  condi¬ 
tion.  Thb  leaves  us  with  two  problems:  (1)  find¬ 
ing  which  elements  to  insert  into  J;  (2)  inserting 
them  with  their  current  weights. 

Suppose  a  maximal  point  b  deleted,  and  let 
(^Li  yi)  and  (xi{,  y^)  be  its  left  and  right  neigh¬ 
bors  in  I,  respectively.  We  introduce  for  con¬ 
venience  virtual  end  elements  for  /,  p.oo  — 
(-oo,-t-oo)  and  p+oo  (+oo,-oo),  so  that  these 


neighbors  are  always  defined.  Any  new  point 
(x,  y)  to  be  added  to  I  satisfies  xi  <  x  <  xr  and 
Vl  >  y  >  Vr-  To  assist  in  finding  all  such  points, 
we  will  have  initially  ordered  the  whole  of  P  in 
lexicographically  increasing  order,  first  on  x,  then 
on  y,  amd  finally  on  the  initial  weights.  Starting  at 
{xR,yR),  we  search  to  the  left  in  the  lexicographic 
ordering  for  the  first  undeleted  point  (ari,yi)  with 
yi  >  Vr-  Starting  at  (xi,yi),  we  search  again  to 
the  left  for  the  first  point  {x2,y2)  with  y2  >  yi. 
We  continue  in  thb  way  until  we  reach  the  point 
(xL'tyL)-  The  points  (xj, yi),  (x2,y2), . . .  obtained 
are  the  new  maximal  points  to  insert  into  I.  A 
data  structure  consisting  of  a  balanced  binary  tree 
with  the  leaves  so  ordered  and  with  the  maximum 
y  coordinate  for  each  subtree  stored  at  the  cor¬ 
responding  internal  node  supports  thb  search  for 
each  successive  point  efficiently.*  Each  new  point 
b  found  in  time  O(logn),  and  deletion  of  any  ele¬ 
ment  from  thb  set  can  be  done  in  thb  time  also. 
The  set  I  is  initialized  by  performing  the  search 
with  (xi,yt)  =  p_oe  and  (xR,yR)  =  p+oc. 

To  calculate  the  current  weight  for  each  new 
maximal  element,  we  use  two  additional  trees  Tr 
and  Tv  to  record  the  H.UPDATEs  and  V_UPDATEs 
respectively.  The  leaves  of  Tr  are  the  pairs 
(xi,ai}  corresponding  to  each  H_UPDATE,  ordered 
by  the  x,-'s.  Internal  nodes  store  the  sum  of  the 
a-values  of  the  leaves  in  their  subtrees.  Simi¬ 
larly,  Tv  records  the  pairs  (y>,6;)  corresponding 
to  V_UPDATEs.  The  current  weight  w  of  a  point 
p  =  (u,  v)  with  initial  weight  wo  is  given  by 

w  =  «^  +  51  +  12  **>• 

Each  sum  can  be  computed  in  time  proportional 
to  the  height  of  the  tree.  Using  balanced  tree  tech¬ 
niques,  the  height  b  O(logm)  after  m  updates, 
and  the  time  to  maintain  the  trees  after  each  up¬ 
date  b  also  O(logm). 

'Kmith  attributes  the  idea  for  such  a  data  structure  to 

McCreigbt  [1|. 


3  Optimal  Assembly  of  Jumps 

3.1  Short  Jump  to  Long  Jump  Con¬ 
version 

Monotonic  priority  sets  have  application  to  a 
problem  of  assembling  code  for  computers  with 
two  sizes  of  jump  instructions.  Short  jumps  ex¬ 
ecute  quickly,  but  the  distance  they  can  jump, 
their  span,  is  limited,  for  example,  to  about  128 
bytes  in  either  direction.  Long  jumps  on  the  other 
hand  are  to  an  arbitrary  destination,  but  are  more 
costly,  both  in  execution  time  and  program  size, 
since  a  long  jump  instruction  needs  more  address 
bytes.  Choosing  which  kind  of  jump  instruction 
to  use  will  normally  be  done  by  the  assembler. 
The  problem  is  a  nontrivial  one  since  making  a 
given  jump  long  increases  the  span  of  any  jump 
traversing  it.  It  is  possible  to  construct  chains 
of  overlapping  jump  instructions  such  that  each 
jump  that  is  made  long  causes  exactly  one  other 
to  exceed  the  maximum  span  for  a  short  jump,  so 
the  latter  too  must  be  made  long. 

It  is  easy  to  see  that  a  minimal  size  program 
can  be  obtained  by  starting  with  all  jumps  short 
and  traversing  the  program  repeatedly,  convert¬ 
ing  short  jumps  to  long  on  each  pass  as  necessary. 
A  straightforward  implementation  of  this  strat¬ 
egy  however  is  quite  inefficient  since  up  to  0(n) 
passes  may  be  required  for  a  program  containing 
n  jumps. 

A  more  conservative  approach  is  to  check  just 
those  short  jumps  within  a  distance  5  of  each  new 
conversion,  where  5  is  the  maximum  span  of  a 
short  jump,  for  these  are  the  only  new  candidates 
for  conversion.  This  approach  leads  to  an  0  (nS) 
time  algorithm. 

Using  a  monotonic  priority  set,  we  obtun  an 
algorithm  for  optimal  assembly  of  jumps  that 
runs  in  time  O(nlogn),  independent  of  5.  A 
jump  from  x  to  y  is  represented  by  the  point 
(-  min(x,y),max(x, y))  with  initial  weight  |a:-y|. 
FIND.MAX  returns  the  remaining  short  jump  of 
largest  span.  It  will  need  conversion  to  a  long 
jump  if  any  jump  does.  Assuming  it  does,  we 
then  delete  it  from  the  priority  set  and  update 
the  priorities  of  the  remaining  elements  by  execut¬ 


ing  V_UPDATE(2,  6)  and  H_UPDATE(-2,6).  Here 
we  assume  that  the  conversion  requires  the  inser¬ 
tion  of  b  extra  bytes  at  position  z.  (Typically, 
z  =  X.)  The  two  given  UPDATE’s  actually  in¬ 
crement  those  jumps  which  span  z  by  26  and  all 
other  jumps  by  6,  preserving  the  priorities  as  re¬ 
quired.  If  we  wish  the  weights  to  equal  the  cur¬ 
rent  jump  spans,  then  a  further  operation  of  say 
V_UPDATE(— 00,  —  6)  could  be  used  to  restore  the 
proper  weights.  The  time  bounds  follow  since  in 
this  case  m  —  0  (n). 

3.2  Multilevel  Jumps 

In  a  logical  extension  to  the  above  problem  we  may 
have  further  levels  of  jumps,  for  example  short, 
long,  and  huge.  This  problem  is  handled  using 
two  monotonic  priority  sets,  one  to  deal  with  con¬ 
versions  from  short  to  long  and  the  other  for  con¬ 
versions  from  long  to  huge.  All  jumps  start  as 
members  of  both  sets.  A  jump  is  deleted  from  the 
first  set  when  it  is  promoted  from  short  to  long, 
and  it  is  deleted  from  the  second  set  when  it  is 
made  huge,  but  no  jump  is  deleted  from  the  sec¬ 
ond  set  until  after  it  has  been  deleted  from  the 
first.  In  slightly  more  detail,  at  each  stage  an  at¬ 
tempt  is  made  to  convert  a  short  jump  to  long, 
and  only  if  that  fails  is  an  attempt  made  to  con¬ 
vert  a  long  jump  to  huge.  All  updates  are  applied 
to  both  sets.  Any  finite  number  of  jump  levels 
can  be  handled  in  the  same  way,  as  can  the  case  of 
short  and  long  jumps  in  which  the  maximum  span 
of  a  short  jump  depends  on  its  direction,  whether 
forward  or  backward  in  memory.  A  continuous 
analogue,  in  which  a  jump  of  distance  d  has  size 
logd,  appears  to  be  much  more  difficult,  and  we 
leave  that  problem  open. 

3.3  >/.r-Completenes8  for  Shared  Jumps 

A  reduction  of  program  length  may  be  possible  if 
we  allow  a  jump  to  a  distant  destination  to  be  as¬ 
sembled  as  a  chain  of  two  jumps,  a  short  followed 
by  a  long,  for  then  several  nearby  short  jumps 
could  share  the  same  long  jump. 

Though  related  to  the  problem  of  Section  3.1, 
optimization  with  shared  jumps  appears  to  be 
very  hard.  Some  microprocessor  codes  have  only 


short  conditional  jumps,  so  that  any  long  condi¬ 
tional  jump  always  requires  two  instructions.  We 
can  show  in  this  case  that  given  a  program  and 
an  integer  k,  the  problem  of  testing  whether  the 
jumps  can  be  assembled  to  yield  a  program  of 
length  at  most  k  is  A//P-complete.  A  proof  of  this 
result  will  appear  in  the  final  version  of  this  paper. 

4  Double  Priority  Queues 

Suppose  we  are  processing  a  set  of  elements  ranked 
by  two  criteria,  for  example,  expected  run  time 
and  memory  size  of  computer  jobs.  The  next  ele¬ 
ment  to  be  processed  and  deleted  is  one  for  which 
the  sum  of  its  two  ranks  among  the  remaining  el¬ 
ements  is  minimal.  This  queueing  procedure  can 
be  implemented  efficiently  using  priority  sets.  An 
element  e  is  identified  with  the  point  (— x,  — y)  in 
a  priority  set,  where  x  and  y  are  its  initial  ranks, 
and  its  initial  weight  is  —x  —  y.  When  e  is  deleted, 
H_UPDATE(-x,-l)  and  V_UPDATE(-y, -1)  are 
performed  to  make  the  weights  reflect  the  sum  of 
the  current  ranks.  Even  though  the  increments 
are  negative,  we  can  prove  that  the  updates  pre¬ 
serve  monotonicity,  so  the  priority  set  implemen¬ 
tation  of  Section  2  can  be  applied. 

5  Conclusion  and  Open  Prob¬ 
lems 

Both  the  monotonicity  requirement  and  the  re¬ 
striction  to  the  plane  appear  essential  to  our  ap¬ 
proach.  Relaxation  of  either  constraint  would 
seem  to  require  new  ideas. 

The  “log”  factor  in  the  time  bound  seems  dif¬ 
ficult  to  eliminate.  For  example,  even  when 
the  points  are  given  in  order  of  decreasing  ini¬ 
tial  weights,  at  least  n(nlogn)  comparisons  may 
be  needed  in  the  case  that  the  weights  can 
be  arbitrary  real  numbers.  To  see  this,  take 
n  points,  pi  =  (l,6i  -  2),...,p„  =  (n,6„  - 
2n),  in  this  order,  where  0  <  6,-  <  1  and 
u;(p,)  =s  6f  -  i,  1  <  «  <  n.  Then  perform 
H_liPDATE(t,  1)  for  each  i  €  {1, . . . ,n}.  Now  the 
resulting  weights  are  and  a  sequence 

of  n  DELETE(F1ND_MAX)  operations  removes  the 


points  from  the  priority  set  in  order  of  decreasing 
value  of  hi,  thereby  effecting  a  sort  of  {6i , . . . ,  6„}. 
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